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THE  EFFECT  OP  TEMPERATUfS  STRATIFICATION 
OF  A  MEDIUM,  ON  TURBULENCE 

A.  S.  Monin 

Institute  of  Atmospheric  Physics,  USSR 
Academy  of  Sciences,  Moscow,  USSR 

In  a  medium  with  turbulent  pulsations  of  density  p',  being  in 

a  field  of  gravity,  a  significant  rol-.^  in  the  dynamics  of  turbulence 

is  acquired  by  Archimedean  forces;  the  accelerations  created  by  them. 

are  directed  vertically  and  have  the  value  -  « where  p  -  averaged 

density  of  medium,  g  —  acceleration  due  to  gravity,  and  the  prime 

signifies  pulsation,  l.e.,  deviation  from  mean  value.  In  the  case  of 

atmospheric  air,  the  consideration  of  which  we  will  limit  ourselves 

to  here,  pulsations  of  density  are  created  mainly  by  pulsations  of 

temperature  and  have  the  form  p'  =  _i7-’,  so  that  the  value  of  Archl- 

^  r 

medean  accelerations  turns  out  to  equal  gy,  where  T  -  average  tem¬ 
perature  (in  sea  water  the  matter  is  more  complicated,  since  pulsation 
of  density  here  Includes  the  contribution  of  pulsation  salinity). 

During  stable  (subadlabatlc)  temperature  stratification  vertical 
shifts  of  air  particles,  created  by  turbulence  of  a  dynamic  origin, 
are  accompanied  by  expenditures  of  energy  for  work  against  Archimedean 
forces  and  lead  to  trans'onnatlon  of  part  of  the  jclnetlc  energy  of 
turbulence  into  potential  energy  of  density  stratification;  conse¬ 
quently,  stable  stratification  weakens  turbulence. 
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During  unstable  (superadlabatlc)  temperature  stratification 
Archimedean  forces,  conversely,  create  vertical  displacements  of  air 
particles  (at  p'  <  0  or  T'  >  0  —  upwards,  at  p'  >  0  or  T'  <  0  —  down¬ 
wards),  which  leads  to  an  Increase  of  kinetic  energy  of  turbulence 
at  the  expense  of  potential  energy  of  stratification  of  density;  con¬ 
sequently,  unstable  stratification  strengthens  turbulence. 

For  description,  of  turbulence  In  a  temperature-stratified  medium 
a  similarity  theory  was  proposed  which  was  founded  on  the  assumption 
that  In  a  layer  with  constant  stress  of  friction  for  height  t  «  pu^ 
(where  u,  —  so-called  friction  velocity)  and  constant  vertical  tur¬ 
bulent  heat  flux  q,  what  Is  the  surface  layer  of  air  with  a  thickness 
of  several  dozen  meters,  all  the  characteristics  of  turbulence  which 
are  not  too  small-scaled  (do  not  depend  on  molecular  viscosity  and 
thermal  conductivity  of  air)  can  depend  only  on  three  constant  measured 
parameters  -  friction  velocity  u,,  parameter  of  buoyancy  g/T,  and 
"temperature  flux"  qlcf\\.  A  systematic  account  of  this  similarity 
theory  is  given  In  Chapter  IV  of  the  book  [1],  where  an  extensive 
bibliography.  Is  also  given. 

In  this  report  we  will  dwell  on  several  separate  problems 
connected  with  the  Influence  of  temperature  stratification  of  the 
medium  on  turbulence. 

Influence  of  Stratification  on  the  Value  of  Vertical 
Turbulent  Fluxes  on  Momentum,  Heat,  and  Moisture 

It  Is  known  that  during  unstable  stratification  and  strong  wind 
turbulence  Is  much  more  intense  than  during  stable  stratification  and 
weak  wind.  We  will  illustrate  this  qualitative  assertion  wlt.^.  some 
quantitative  findings.  We  will  examine  a  situation.  In  which  the 
height  of  roughness  of  the  underlying  surface  Zq  Is  considered  known 
and  gradient  measurements  are  made  on  three  heights  z  »  H/2,  H,  and 
2H.  Measured  primarily  here  are  the  average  wind  velocity  u  u(H) , 
difference  of  temperatures  «T  -  T(2H)  -  T(H/2),  and  difference  cf 
values  of  specific  humidity  -  Q(2H)  -  Q(H/2).  Then  friction 
velocity  u,,  vertical  turbulent  heat  flux  q,  vertical  turbulent 
moisture  flux  E  (rate  of  evaporation)  according  to  the  above-mentioned 
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aimllarlty  theory  can  be  presented  In  the  form 


(1) 

(2) 

(3) 


where  R  —  Richardson  satirical  number,  determined  by  the  formula 


(H) 


Further  for  concreteness  we  will  consider  that  H  ■  1  n  and  Zq  •  1  cm. 
In  cases  of  8tz>ong  Instability  R  has  a  value  of  the  order  -0.05,  and 
In  cases  of  strong  stability  -  of  the  order  to.  3  (during  neutral 
stratification,  of  course,  R  ■  0).  Theoretical  and  empirical  data 
concerning  functions  f,  4,  and  ip  (see  [1]  and  [2])  show  that  function 
f,  characterizing  the  ratio  u,/u,  during  strong  Instability  has  a 
value  of  around  0.11-0.12,  at  neutral  stratification  around  0.09,  and 
at  strong  stability  around  0.06-0. 07,  so  that  during  Instability  It 
turns  out  to  be  one  and  a  half  or  tyo  times  greater  than  during 
stability.  Functions  t  and  i|i  change  considerably  more  sharply  -  with 
a  change  of  R  from  -0.05  to  +0.3  they  decrease  by  25  times. 

Values  iqi  and  |E|  change  still  more  sharply  due  to  the  fact 
that  between  u  and  {T  (and  apparently  also  between  u  and  5Q)  there 
la  a  correlation;  strong  temperature  Inversions  (positive  6T)  are 
observed,  as  a  rule,  with  a  weak  wind,  but  during  strong  Instability 
(negative  iT),  conversely,  there  Is  almost  no  calm.  As  an  example 
the  figure  represents  a  graph  of  correlation  between  u  and  6T  for  the 
period  July-September  1959  In  the  area  of  Tslmlyansk  according  to 
work  [33,  which  Indicates  that  values  for  the  product  of  lu6r|  during 
Instability  usually  several  tlm-ss  greater  than  during  stability 
(typical  example  of  Instability  —  u  «  6  m/s  and  4T  •  -1.5®;  typical 
example  of  strong  stability  -  u  »  0.5  m/s  and  ST  -  3°;  the  product 
of  |u4T|  here  changes  by  six  times).  As  a  result  values  IqI  and  |E| 
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during  strong  Instability  turn  out  to  be  two  orders  greater  than  dur- 

Ing  strong  stability,  so  that  the  influence  of  strati - - - 

turbulent  heat  transfer  and  passive  impurities  turns  out  to  be  m.,nh 
stronger  than  its  influence  on  turbulent  transfer  of  momentum. 


The  strong  dependence  of  values  f,  ♦,  and  on  the  parameter  of 
stratification  R  and  the  presence  of  a  correlation  between  u  and  6T 
(and,  apparently,  between  u  and  iQ)  prevent  the  calculation  of 
average  climatic  values  of  turbulent  fluxes  of  momentum,  heat  and 
moisture  based  on  average  climatic  values  u.  «T,  and  jQ  with  the  help 
of  formulas  of  the  type  (l)-(3);  such  a  calculation  could  lead  in 
certain  cases  to  errors  of  an  o^er  of  hundreds  of  percents  [3]. 


Average  quadratic  pulsations  of  velocity  have  values  of  the 
order  u,  -  uf  and,  consequently,  during  instability  will  usually 
be  only  several  times  larger  than  during  stability.  Pulsations  of 
temperature  have  values  of  the  order  r.-r.lV.srA.  and  during 
Instability  they  will  be  an  order  greater  than  .during  stability. 

Let  us  note  that  significant  pulsations  of  temperature  can  appear 
inside  a  turbulent  layer  only  due  to  vertical  mixing  at  averaged 
potential  temperature  e(s)  which  is  changing  with  height.  At  e(s)  - 
-  const,  l.e.,  in  the  case  of  neutral  stratification,  pulsations  of 
temperature  will  not  appear,  they  can  only  get  into  the  examined 
layer  from  without,  for  example,  in  the  surface  layer  of  air  -  from 
thermal  heterogeneities  of  the  underlying  suface.  But  the  latter 
change  in  the  course  of  twenty-four  hours;  by  day  the  dark  spots  are 
wanner,  and  at  night  this  is  not  so.  It  is  possible  to  think  that 
neutral  stratification  corresponds  exactly  to  periods  of  levelling 
off  of  thermal  heterogeneities  of  the  underlying  surface.  They  during 
neutral  stratification  in  general  there  will  be  no  significant  pul¬ 
sations  of  temperature  and,  consequently,  a  situation  appears  which 
at  the  first  sight  seems  paradoxical:  during  transition  fron 
to  stable  stra(;lfication  pulsations  of  temneref.n^  ...ould  mcree,, . 
^thqugh  turbulence  (i.e.,  pulsations  of  velocity)  le  suporeseed 
These  measurements  of  pulsations  of  temperature  during  various  thermal 
stratification  show  [*.]  that  this  curious  situation  Indeed  takes  place. 
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u 
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Correlation  betaeen  u  and  <T  for  the 
Septenber  1959. 


'  S>*T 

period  June- 


Just  aa  in  the  case  of  a  teaperature-unlform  tedium  it  la  natural 
to  assume  that  turbulence  In  a  temperature-stratified  medium  will  be 
jpally  un^fote.  This  means  that  conditional  probability  distributions 
lor  differences  of  ealues  of  hydrodynamic  fields  in  sufficiently  close 
points  of  space  of  time,  under  the  condition  that  rate  of  movement 
Of  the  origin  of  coordinates  0  of  the  chosen  Inertial  system  of  ref¬ 
erence  S  Coincides  with  Initial  Euler  speed  u.  .  u(0.  t„)  at  point  0 
can  depend  only  on  space-time  coordinates  of  points  of  observation  In 
system  S,\but  do  not  depend  clearly  on  where  x„  ('-.i.  - 

coordln^tjjfs  of  point  0  at  initial  moment  of  time  t^. 

In  particular,  distributions  of  probabilities  fcr  differences 
of  meteorological  elements  In  close  points  can  clearly  depend  only  on 
differences  of  heights  of  point,  of  observation,  and  on  full  height 
z  of  one  of  them  can  depend  only  Impilclty  (according  to  the  hypotheses 
of  slml.’arlty  of  A.  N.  Kolmogorov  and  A.  M.  Obukhov  -  only  by  moans 
of  average  velocity,  dependir«  on  z,  of  dissipation  of  turbulent 
energy  c  and  average  rate  of  balancing  of  temperature  heterogeneities 
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Inequality  of  different  directions  In  space,  created  by  the  pre¬ 
sence  Of  average  flow  and  limiting  fluxes  of  walls  manifest  them¬ 
selves  on  the  properties  of  large-scale  components  of  turbulence-,  but 
8houlc<  not  Influence  the  properties  of  its  small-scale  components  and 
the  local  characteristloB  of  turbulence  determined  by  them.  However, 
in  a  tenperature-stratlfled  mcdl urn  all  fluctuations  of  density,  both 
large-scale  and  small-scale,  will  experience  the  influence  of 
Archimedean  forces.  Therefore  the  vertical  direction  will  be  assigned 
for  components  of  turbulence  of  any  scales,  and,  consequently,  tur¬ 
bulence  In  a  temperature-stratified  medium  cannot  be  locally  Isotropic 
[5].  But  it  is  natural  to  assume  [6,  7J  that  It  will  be  locally 
axlaymmebrlcal  with  respect  to  vertical,  i.e.,  that  the  probability 
dlKtrlbt.tlon  for  differences  will  be  invariant  relative  to  rotations 
of  the  system  of  coordinates  around  a  vertical  axis  and  specular  re¬ 
flections  In  any  vertical  planes. 

3-  Similarity  Hypotheses  For  Inertlallv  Convective 
Interval  of  Turbulence  wectrum  in  a 
^Stably  Stratified  Medlufi 

According  to  the  known  hypotheses  of  similarity  of  A.  K.  Kolmogorov 
and  A.  M.  Obukhov,  the  statistical  characteristics  of  components  of 
turbulence  with  scales  from  the  so-called  inertial  Interval  (small  in 
comparison  to  external,  but  large  In  comparison  to  internal  scale  of 
turbulence)  in  a  case  of  neutral  stratification  of  the  medium  can  de¬ 
pend  only  on  the  two  measured  parameters  e  and  N  which  were  mentioned 
above.  According  to  A.  M.  Obukhov  [8],  in  the  case  of  a  temperature- 
stratified  medium  one  should  add  to  them  the  parameter  of  buoyancy  g/T 
which  makes  it  possible  to  introduce  a  new  scale  of  length. 


If  this  scale  belongs  to  the  inertial  Interval,  then  in  such  a  case 
the  Interval  of  scales,  small  In  comparison  to  external,  but  large  in 
comparison  to  internal  scale  of  turbulence,  can  naturally  no  longer 
be  called  simply  Inertial,  but  an  inertial-convection  interval. 
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In  vase  of  veaiy  strong  atSbility  a  large  share  of  energy  of 
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'■fdssllte  to,  'psaastt  of.-  icaases  to .  be 

■alfn^ftcant^for  .ebni^^lonii  ’4i‘''turtttieiice''l,ij,,  >Sbe  oOnveptioa -'Interval:. - 
(oonatltutlhi'  or  ^l»-'l^rUal,con»eetlob ' 

Interval),  so  that  .Migl^ri  here  li^.de.tensliied  .ohjly  by  parfiaetaisi 
I  and  Then  spS^M'density  of  ene:^  i[(K)  and  apectral  density  ; 


of,  heter< 
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(7) 

(Bolglano  [43,  «o»iinrt^P  another 

hypothesis  for  de8crlbl%il^  eoat^idn  intSrval  of  the  spectrum  of 
turbulence  in  the  case  of  stro^  >«iAJXU#,  '--  aR'.sh  i»* 

reglmea  of  turbulende  in  Jthls  Interval  ds  determined  by  the  paiweter 
of  buoyancy  g/T  Snd  grableNdt  of  aveMge  potential  tenipeip'ature  ""r.  so  that 
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If  we  tabs  the  valuessf -g/T  and  r  for- charaCterlstir'dimensibn  '  par¬ 
ameters  also  for  the  spectrum  of  temperature,  then  frbm  consiieratioM 
of  dimension  we  .obtain-  ,-  ^  ' 
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The  etrongest  differenci  here  Is  revealed  In  the  spectra  of 
temperature  field  (7)  and  ('»);  unfortunately,  so  far  we  do  not  have  t 
experimental  data  which  woiild  make  It  possible  to  select  one  of  these 
spectra  and  reject  the  other.  Let  us  note  only  that  spectrum  (9) 
turns  out  not  to  be  dependent  on  the  parameter  of  buoyancy,  which 
seems  strange.  Let  us  stless,  mhChefeoro ,  that  thm  assumption  ccn^- 
cemlng  the  dependence  of  conditions  of  turbulehcs  In  the  convection 
Interval  of  the  spectrum  on  the  gradient  of  >'ave‘F*£e  "'temperature,  but 
not  on  the  gradient  of  average  veloc.lty  (enteri|^  into  the  expression 
of  the  work  created  by  Reynolds  stresses) f  senui’  liioglcal  -  it  Is 
natural  to  assume  that  parametaia  and  f.  *  ib/.v  should  enter  symmet¬ 
rically  Into  tlhe  expressions  for  statistical  characteristics  of  tur¬ 
bulence. 

f-'  ‘  ... 

■  tlmttlMC  •Instability  and  Limiting  Stability 


In  description  of  turbulence  In  llmltlnjg  eases  of  strong  Insta¬ 
bility  and  \^ikt^«auBrsp  a*^  Let 

start  with  the  case  of  strong  Instability.  According  to  the  previously 
^***.®*T  ^  tbf.  Sharacteplstlcs  of  turbulence  In 

ito  dspop4  00, ratf;  of  motion  u,and  mre  dstomlned 
,pnly  .by-|^s  two  pajpamstefs  gA  sn4  e/op.  so  that;  .for.axsapie,  the 
,  gradient  of  average  potential  temperature;  a  in  the 'Surfaoe:»  layer  of 
alf  turns  out  to  be  dependent  pn  height  according  to  the,  law 


f .  r\ 


(10) 


(although  ao/Oz-  O  ^jth  Increase  of  zi  hire  It  Is  quite  Impossible  to 
consider  that  temperature  etrmtlflcatlon  of  medium  approaches  neutral, 
since  the  Richardson  number  with  an  Ineivase  r>f  fe  aspires  to  and 
never  to. sero).  On  the  other  hand,  NalkuS  [133  constructed  theory  of 
thermal  convection,  accord.ing  to  which 'flo/ai  - i-*  irith  a  proportionality 
factor  depending  on  the  coefficient  of  moleieUlar  thermal  conductivity  ' 
of  air  (the  latter  ls  somewhat  strange  for  turbulence  developed  during 
conveotloni  a  hypothetical  explanation,  proposed  for  cases  With  the 
absence  of  wind  and  created  mainly  by  wind  from  horizontal  mixing,  * 
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Involves  the  influence  on  ascending  convection  streaos  of  conditions 
of  their  conception  in  th',  underlayer  of  molecular  thermal  conduction) . 

Empirical  data  convlndlhgly  show  that  formula  (10)  during  un¬ 
stable  stratification  is  Justified  at  least  in  the  layer,  in  which 
0.02  <  |e/L|  <  1  (where  L  -  known  scale  of  length,  composed  of 
parameters  g/T,  and  u,).  However,  there  are  some  measurements 

taken  in  the  atmosj^hbre  (Webb)  and  under  laboratory  conditions 
(Thomas  and  Townsend),  according  to  which  at  |z/L|  >  1  formula  (10)  no 
longer  is  conflMed,  and  perhaps  even  the  conditions  are 

established.  A  dlaeosslon  ef  this  difference  can  be  found  in  section 
8.2  of  book  [1],  where  the  corresponding  literature  is  cited.  Let  us 
point  out,  furtheiaere,  the  review  article  by  Spiegel  £13]  and  the 
more  recent  Woi4t  df  Dyer  [l<t],  in  shlch  on  the  basis  of  the  latest, 
very  thorough  measurements  in  Kerang  (North  Victoria,  Australia)  a 
definite  conclusloh  fs  made  in  fivof  of  foiwula  (itf)  (yjth  at  |z/L|  <  l 
and  at  |z/Lt  >1.  j  * 

Still  more  difficult  is  the  case  of  limiting  stability.  The 
differences  appearing  here  are  formulated  most  conveniently  in  terms 
of  the  ratio  a  ■  K^/K  of  coefficients  of  exchange  for  heat  and  for 
momentum.  This  value,  which  at  neutral  stratification  is  arcond  a 
unit,  apparently  decreases  with  em  Increase  of  stability,  and  the 
question  arises,  does  it  strive  within  a  limit  to  a  certain  positive 
limit  cigp  >  0  or  to  zero.  The  so-called  dynamic  Richardson  number 
Rf  -  oRl  (  where  Hi  -  ordinary  Richardson  number)  with  an  Increase 
of  stability  aspires  to  a  finite  limit  Rf  .  Consequently,  in  the 
case  “  “jp  >  0  the  ordinary  Richardson  number  will  aspire  to  finite 
limit  a'ld  in  the  case  a  -►  0  will  Increase  without  limit.  In 

the  first  of  these  cases  according  to  the  similarity  theory  ([1], 
point  7.3)  the  gradient  of  temperature  <>T/di  will  asymptotically 
approach  a  constant 


^ _ »  »  / 


«,p».  1 . 
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BO  that  temperature  will  Increase  linearly  with  height:  in  the  second 
case  tenperature  will  Increase  with  height  more  rapidly  than  by 
linear  law. 

Speaking  In  favor  of  the  assumption  a  *  a  >0  are  the  observa¬ 
tions  of  Llljequlst  [15],  who  detected  linear  profiles  of  temperature, 
corresponding  to  formula  (11),  during  strong  inversions  In  Antarctica, 
and  the  results  of  processing  the  empirical  findings  of  various  authors 
by  A.  B.  Kazanskiy  [16],  according  to  whom  R1  <  0.3.  Let  us  mention 
also  our  works  [7]  and  [1?],  In  wWcA^wlth  the^help  of  simplified 
Prledmar-Keller  equations  for  sl^ie-polnt  second  moments  of  velocity 
of  Hn4  and  teiAperature  we  obtalrid  the  depen<icnoe  of  a  on  Hf,  leading 
to  the  limiting  value  a  >  0  (In  an  example  from  [17]  a  \  0.2  and 

HI  , Ik  0.3).  .  up  .  .  .  Kp  ^ 

KP  ' 

On  the  other  hand,  there  are  certain  measurements  In  nature 
(J.  Ti^rlor-Prowdmenj  Koleanlkov)  and  in  the  laboratory  (Ellison  and 
Turner)*  a  discussion  of  can  bi  foun^i;  In  section  8.2  book  ■ 

[1],  which  are  In  qualitative  agreement  with  the  formula  of  Ellison 

f  —  R]/Rf  '■ 

kHen  (12) 

(the  oor.ciuslph  of  which  la  bas4a,  it  la  true,  on  much  rougher 
almpllflcatlons  of  Frlednan-keller  equations  than  those  which  are 
used  In  [7]  and  [17]).  Based  on  these  findings  Rf  j,  0.10-0.15  and 
o  %  0.02-0,05  at  A  <  Rl  <  10. 


For  the  solution  of  this  difference  in  the  Values  of  a  and 
f“*^her  accumulation  of  findings  concerning  turbulence  under 
conditions  of  very  strong  stability  Is  necessary.  Besides  one  should 
consider  that  under  such  conditions  turbulence  Is  apparently  concen¬ 
trated  in  areas  which  are  randomly  distributed  In  space  and  with 
vertical  dimensions  which  are  much  less  than  horizontal.  Inside 
such  regions  turbulent  mixing  leads  to  vertical  temperature  balance 
(more  exactly  -  potential  temperatures) ,  and  outside  their  temperature 
can  change  very  sharply  with  height.  As  a  result  Instantaneous 
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profiles  of  temperature  can  have  numerous  randomly  distributed  breaks, 
and  only  a  sufficiently  wide  averaging  will  lead  to  smooth  average 
profiles  of  temperature,  perhaps  of  the  form  (11).  Under  such 
conditions  It 'Is  expedient  not  to  be  limited  to  a  determination  of 
only  average  profiles,  insufficiently  characterizing  the  structpre  of 
turbulence,  but,  by  following  methods  of  describing  alternating  tur¬ 
bulence  for  the  borders  of  turbulent  boundary  .layers,  trails,  and 
streams,  to  determine,  for  example,  the  field  of  "coefficient  of 
intermit tenc'e,"  l.e.,  the  probability  of  the  presence  of  tuirbulence 
In  different  points  of  space. 
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abstract 


(U)  following  problems  associated  with  the  effect  of  tem¬ 
perature  stratificatlcKi  on  turbulence  are  discussed;  1)  the 
effect  of  stratification  on  the  vertical  turbulent  fluxes  of 
momentum,  heat  and  moisture;  2)  the  symmetric  properties  of  turbu¬ 
lence  in  a  temperature  stratified  medium;  5)  similarity  hypotheses 
for  the  inertlally  convective  interval  of  turbulence  spectrum 
in  a  stability  stratified  medlvun;  4)  the  limiting  cases  of 
instability  and  stability.  The  influence  of  temperature  stratifi¬ 
cation  and  turbulence  is  described  with  the  aid  of  three  parameters- 
the  buoyancy  parameter,  the  friction  velocity  and  the  so-called 
'temperature  f lux" . (  JOrig.  art.  has;  1  figure,  12  formulas. 


